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Abstract: 

We employ nonperturbative flow equations for the description of the effective action 
in Yang-Mills theories. We find that the perturbative vacuum with vanishing gauge 
field strength does not correspond to the minimum of the Euclidean effective action. 
The true ground state is characterized by a nonvanishing gluon condensate. 



1 Introduction 

A perturbative calculation of the effective action for constant colour-magnetic fields 
in Yang-Mills theories indicates that the configuration of lowest Euclidean action 
does not correspond to vanishing fields |]ll)0- This has led to many interesting 
speculations about the nature of the QCD vacuum. Unfortunately, perturbation 
theory is clearly invalid in the interesting region in field space. It breaks down 
both for vanishing magnetic fields and for fields B corresponding to the minimum of 
the perturbatively calculated effective action T[B]. Despite many nonperturbative 
indications for a nontrivial QCD vacuum from lattice studies and phenomenology, 
an analytical establishment of the phenomenon of gluon condensation is still lacking 
so far. 

In this paper we present new analytical evidence for gluon condensation based 
on the non-perturbative method of the average action [^. The average action Tk 
is the effective Euclidean action for averages of fields which obtains by integrating 
out all quantum fluctuations with (generalized) momenta g^ > k"^. It can be viewed 
as the standard effective action F computed with an additional infrared cutoff ~ k 
for all fluctuations. In the limit k ^ the average action equals the usual effective 
action, Fq = F. For A; > no infrared divergences can appear in the computations 
and a lowering of k allows to explore the long-distance physics step by step. The k- 
dependence of the average action is described by an exact nonperturbative evolution 
equation [Q.^ The structure of this equation is close to a perturbative one-loop 
equation but it involves the full propagator and vertices instead of the classical ones. 
For gauge theories it can be formulated in a way such that F^ [A] is a gauge-invariant 
functional of the gauge field A |0].Q 

The flow equation for the pure non-abelian Yang-Mills theory has been solved |^ 
with a very simple approximation - the average action F^ has been truncated with a 
minimal kinetic term ~ ZkF^^F^^^ . From the /c-dependence of Z^ the running of the 
renormalized gauge coupling g{k) has been derived for arbitrary dimension d. Most 
strikingly, this lowest-order estimate suggests that the wave function renormaliza- 
tion Zfc reaches zero for k = k^o and turns negative for k < k^o- Here the scale koo 
can be identified with the confinement scale, i.e. the scale where the renormalized 
gauge coupling g'^{k) ~ Z^T^ diverges. This was a first nonperturbative analytical 
indication for the instability of the perturbative vacuum with F^^j = 0. Indeed, 
negative Z^ implies gluon condensation - the minimum of F^ must occur for a non- 
vanishing gauge field for all fc < fcoo Pi- It is obvious, however, that an effective 
action ~ ZkB^ is insufficient to describe gluon condensation. One expects positive 
Ffc for large B and this cannot be accomodated with a negative Zk in the "lowest 
order truncation" . In this paper we enlarge the "space of actions" by considering for 
Ffc an arbitrary function of constant magnetic fields, i.e. F^ ~ Wfc(|i?^). The flow 
equation describes how the function Wk changes its shape, starting from a linear 



"'^For the relation to earlier versions of exact renormalization group equations H see ref. [0. 
^See ref. 01 for alternative formulations. 



dependence Wk = \ZkB'^ for large k. We will see that within our truncation Wk 
develops its absolute minimum for 5 7^ once k becomes smaller than a critical 
scale kci and that the linear term indeed vanishes for koo < kc- Despite some short- 
comings of the present truncation for k^k^o (which are discussed in the conclusions) 
we interpret these findings as a strong analytical indication for gluon condensation! 
In the following we shall consider the pure SU{N) Yang-Mills theory. For this 
case the exact evolution equation of ref. reads (t = In k) 
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{t^^\a, A] + Rkiv^^^A, A])y ^i?.(rf [AA]) 

-D>^[A]D,[A] + Rk{-D'[A]))-'^^Rk{-D'[A]) 
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Here the first trace on the r.h.s. arises from the fluctuations of the gauge field. It 
involves on integration over space-time ( "x" ) as well as a summation over color ( "c" ) 
and Lorentz ("i^") indices. The "color" trace is in the adjoint representation. The 
second trace is due to the Faddeev-Popov ghosts and is over space-time and color 
indices only. The precise form of the infrared cutoff is described by the function Rk- 
It is convenient to choose 



Rk{u) = u[exp{Zf^ u/k ) — 1] 
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In general the wave function renormalization Zk could be a matrix in the space of 
fields which may even depend on A. (In ref. we used Zk = 1 for the ghosts and 
a /c-dependent constant Zk = Zp^k for all modes of the gauge field.) Furthermore, 

/r)\ — 

r^ [A, y4] denotes the matrix of second functional derivatives of F^ with respect 
to A, with the background field A kept fixed. The modes of the gauge field are 
declared "high-frequency modes" or "low-frequency modes" depending on whether 
their eigenvalues with respect to the operator F^ '[A, A\ = F^ '[A, A\ |^^^ are larger 
or smaller than /c^, respectively. Therefore it is this operator which appears in 
the argument of Rj^ and in this sense Rk acts as an effective infrared cutoff by 
suppressing the low frequency modes. The flow equation ( |1.1| ) can be rewritten in 
close analogy to a one-loop formula 
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■^In H we used the classical 5^^-' [A, A\ rather than T\. [A, A\ for this purpose. While preserving 
all the general properties of Tk, the new flow equation is much easier t o ha ndle from a technical 



point of view. This change also entails the different positioning of Zk in (L2) relative to the one in 
0: for the simple trun cati on used there one has F^ ' = Zp.kS'-^^ for part of the modes. We only 
will use the definition (1.2) if for high momenta g'^ — > 00 one has u — > 00. We should mention, 
however, that for negative eigenvalues of Fj. [A, A] the vanishing of Rk for A: ^ is guaranteed only 
if the ratio u/k"^ in (1.2) remains finite. This problem concerns mainly the approach to convexity 
of the effective action for /c — > and is of no relevance for the present work. 



The derivative -^ acts only on the exphcit ^-dependence of the function Rk, but not 

on r^ [A, A] . It is now easy to describe the relation between the effective average 
action Tk and the conventional effective action. Let us first make the approximation 
-^ -^ ^ in eq. ( |1.3| ). This amounts to neglecting the running of Tk on the r.h.s. of 
the evolution equation. It is then trivial to solve it explicitly: 
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-D^[A]D,[A] + Ra{-D'[A])] } + O (^|rf ) 
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with A some appropriate high momentum scale (ultraviolet cutoff) where we may 
identify Fa with the classical action S including a gauge fixing term. This formula 
has a similar structure as a regularized expression for the conventional one-loop 
effective action in the background gauge [|13, 0]- There are two important differences, 
however: (i) The second variation of the classical action, S^'^\ is replaced by F^ . 
This implements a kind of "renormalization group improvement" and transforms 



( |1.4| ) into a sort of "gap equation", 
explicit infrared cutoff Rk- Because 



(ii) The effective average action contains an 



lim Rkiu) 
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a /c-dependent mass-type term is added to the inverse propagator F^ for low fre- 
quency modes (m — > 0), but it is absent for high frequency modes {u — > oo). For the 
purpose of comparison we will also consider in this paper a second choice for the 
cutoff, namely simply a constant 



Rk 



Z.k'^ 
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Despite the similarity of ( |1.4|) with a one-loop expression, we stress that, for 
k —>■ 0, the solution of the original renormalization group equation ( |1.1| ) or (|1.3|) 
equals the exact effective action which includes contributions from all orders of the 
loop expansion. For a detailed discussion of the approximation ( |1.4| ) in the case 
of the abelian Higgs model we refer to [11|, and to ref. 
nonperturbative evolution equations. 



12] for the corresponding 



2 Truncating the space of actions 

In order to find nonperturbative approximative solutions of (|1.1|) we employ the 
following ansatz for F^: 

rk[A,A]=Jd'x Wk{\F^^{x)Fr{x)) + ^Jd''xY.{DMW-An)l (2.1) 



Here Wk is an arbitrary function of the invariant ^F^ with F^^^ the field strength of 
the gauge field A. The fc-dependence of Wk will be determined by inserting (|2.1| ) into 
the evolution equation (|1.1|) . If we think of Wk{0), 6 = jF^^^F^'-', as a power series 
in 6 our ansatz contains invariants of arbitrarily high canonical dimension. Since all 
invariants which occur are of the form 6\ only the dimensions 41, 1 = 1,2, ..., actually 
occur. Also for a fixed dimension 41 the truncation ( ^.11 ) does not contain a complete 
basis of operators. Nevertheless one may hope that ( p.l|) gives a qualitatively correct 
picture of the physics in the regime where the renormalization group evolution has 
already drastically modified the classical Lagrangian ^F"^. We remark that effective 
actions which depend on 6 only play also a central role in the leading-log models 
[|T^ of QCD. The second term on the r.h.s. of ( p.l[ ) is a standard background gauge- 
fixing term [|l^], [^ with a fc-dependent gauge-fixing parameter a^. In addition to 
the ^-dependence of the function Wk{9) we should, in principle, also compute the 
fc-dependence of a^- We will omit this here since the general identities of ref. 
imply that a is independent of fc in a first approximation. 

We should mention at this place that a truncation is actually not defined by 
the terms retained but rather by specifying which invariants in the most general 
form of Ffc are omitted. One may parametrize a general Tk by infinitely many 
couplings multiplying the infinitely many possible invariants which can be formed 
from the gauge fields consistent with the symmetries. In the corrresponding infinite 
dimensional space a truncation is a projection on a subspace which is defined by 
setting all but the specified couplings to zero. (In our case the subspace remains 
infinite dimensional.) In practice, we will choose a particular test configuration A^ 
corresponding to a constant magnetic field. Its detailed definition will be given 
below. The truncation should then be understood in the sense that we use a basis 
for the invariants where all invariants except those used in ( f^.l[ ) vanish for the test 
configuration. By putting the coefficients of all invariants which vanish for the 
test configuration to zero the truncation is uniquely defined. A computation of 
Wk therefore amounts to a computation of the fc-dependent effective action for a 
(particular) constant magnetic field. 

Upon performing the second variation of the ansatz (|2.1|) , 

6'Tk[A,A] = jd''x5A^yVf[A,A]y;jA: (2-2) 

we arrive at the following operator V), : 
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oik 
with 



r?^[^, A]% = W;,{9){Vt[A] - V4A])ll + W;!{9)Sy:[A] + ^iV,[A])ll (2.3) 



= \F^.Fr (2.4) 

Here we used the notation {w, y, z are adjoint group indices and g is the (bare) gauge 
coupling) 



{Vl)% = -{D®D)ll = -Dy-Dr 

S',t = Py^F^^iD^DT (2.5) 

with the covariant derivative {Dy\A\)y^ = 3^5^^ — igA^^{Tz)y^ in the adjoint rep- 
resentation and F^^ = ^^^(Tj)^"'. Moreover, Wj. and Wj! denote the first and the 
second derivative of Wk with respect to 6^. In writing down eq. ( p.4|) we made the ad- 
ditional assumption that the field strength -F^,/[A] is covariantly constant F^^.^ = 
or 

[D,[A],F^,[A]]=0 (2.6) 

In fact, we may choose as a particularly convenient test field a covariantly constant 
color-magnetic field [0 with a vector potential of the form 

A^(x) = n^A,{x) (2.7) 

Here n^ is a constant unit vector in color space (n^n^ = 1), and A^(x) is any 
"abelian" gauge field whose field strength 

^f,u = df^Au - dyk^ = Be^^ = const (2.8) 

corresponds to a constant magnetic field B along the 3-direction, say. (We de- 
fine e^2 — "^21 — I5 with all other components vanishing.) Hence we have 6 = 
\F^jjF^'^ = |i?^ such that 9 and W^. etc. commute with all operators. It is easy to 
see that (p.7|) and (p^) obey the condition (|2.5|). The choice (p.7|) has the advantage 



(2) 
that it allows for an explicit diagonalization of the operator F^ . We finally note 

that Wk can be extracted from rfc[74,y4] which is a gauge-invariant functional of A 

— /o\ 

obtained by putting A = A. It is therefore sufficient to know F^ [A, A]. 

Before returning to the evolution equation, we list a few special properties of the 



covariantly constant fields, which will prove helpful later on. From (p.5|) it follows 
that Afj^ satisfies the classical Yang-Mills equations D^F^y = 0. This in turn is 
sufficient to prove that the operators Vt and Vi commute. As a consequence, one 
may define generalized projection operators [0 

Pt = 1-Pl (2.9) 

which satisfy P^^ = Pt,l, Pt + Pl = I and PtPl = = PlPt- For A^ = they 
reduce to the standard projectors on transverse and longitudinal modes: 

(Pf )^. = d,d,/d' (2.10) 

Furthermore, if A^ is of the form ( p.7| ), it is natural to define another pair of orthog- 
onal projectors, 

pyz ^ ^yz _ ^y^z^ pyz ^ ^y^z^ ^2.11) 



which project on the spaces perpendicular and parallel to n^, respectively. For the 
vector potential (|2.7] ) the matrix A^T^, reads in the adjoint representation 

Ay;{x) ^ {A;T^)y^ = ify-^n^k,{x) (2.12) 

The antisymmetry of the structure constants /^"'^ implies that P\\ and P± commute 
with D^, D^, Vt, Vl and F^^, and that 

P\\A^ = 0, P\\D^ = P\\d^, P\\{VT)^.u = -d^6^.P\\ (2.13) 

The operator S from ( p.5|) factorizes according to 

cyz _ pyz 

'-'nu ~ -'II ■'M'^' 

S,u = ^^.p^a.dPd" (2.14) 

Hence S commutes with D^,Vl and Vt because it annihilates the gauge-field con- 
tained in these operators: 

SD^ = Sd'^, SVT = -Sd'^, SD®D = Sd®d (2.15) 

In physical terms this means that those components of the gauge fluctuations 
SA^^ = a^ which are parallel to n^ decouple from A^ to some extent. In fact, in 
terms of the projections a^'" = P±,\\a^ the quadratic action (|2.2|) with (|2.6D reads 



^fii/\0' 



6'Tk[A, A]= f d'^x {^^^[-d'Wl^S^, + (W;, - —)d^d, + Wis, 

J I OLk 

+ai^mVr + {-- W[)V,]llaj'^] (2.16) 

We observe that the a"-modes couple to the external field only via the derivatives of 
Wk = Wki^B^). In a conventional one-loop calculation one uses the classical Yang- 
Mills Lagrangian \F^j^ rather than Wk{\F^j^). In that case the quadratic action for 
the small fluctuations is given by (|2.16|) with W^. = 1 and Wj! = 0. Hence the one- 
loop determinant resulting from the integration over a" is field- independent and may 
be ignored. In the present case, the a"-modes are important for the "renormalization 
group improvement" , however. The quadratic form (|2.16|) can be diagonalized even 



further by introducing the longitudinal and transversal projections 

all'^ = P,a\\ = Pf all, a^'^ = P^a^ = P^a^ 

a^'^ = P^a^, a^'^ = Pra^ (2.17) 

By virtue of Pf ^s = sPf ^ = 0, P^°^s = sP^ = s and [PLiT), 0] = one obtains 
6^Tk[A,A] = Jd''x{4''''^[-d'w;,6,, + W'ls, 



^fii/\0' 



,Tm,z 



+ ±aH,i,M[_52]^||,L,. 

+ -a^^'^''^[PT]>f '^'^ ^ (2.18) 



(2) 

This block-diagonal form of F^ will facilitate the evolution of the traces occurring 
in the evaluation equation. For example, a"'^ gives no A-dependent contribution 
and, except for an irrelevant constant, the only dependence of F^ on ak arises from 
a-*-'^. Writing 

F(2) [A, A] = fP + Ff + Ff + Ff (2. 19) 

where, in an obvious notation 



(2.20) 





Ff) = P\\PtV'^^IP\\Pt, F^') = P\\PlV^^IP\\Pl 




rf = P^PtT%P^Pt, Ff = P^PlT^I\p^Pl 


with [P\i^,Pl,t\ 


= and 




Vfvf = for A^B 




rp{2) p(2)T _ r. 



(2.21) 



one obtains ([rn ^, -f^i-Pr] = etc.) 

i?fc(F(2)) = P\\PtRu{T^\)P\\Pt + P\\PLRk{Tfj)P\\PL 



+P^PTRkiT%)P^PT + P±PLi?fe(rjyPxPL (2.22) 

3 Evolution equation for Wk 

In this section we derive an evolution equation which governs the scale-dependence 
of the function Wk- We will choose the matrix Z^ in the definition of Rk ( |1.2| ) as 
Zk = 1 for the ghosts and 

Zk = ZkPT[A] + ZkPL[A] (3.1) 

for the gauge boson degrees of freedom. Here Zk, Zk are fc-dependent constants and 
we observe that the choice ( |3.1| ) is compatible with ( [2.221 ). If we insert the truncation 
(13) into (13) with A = A, we obtain {6 = ^B^) 



^§-Wkie) = ^Tr^^,[HiT^\A,A])] 

+ lTT^.L[P±PL{H{Tfh - H{T^ll))] (3.2) 

where F^ is given by (|2.6|) and fi = / d'^x. In ( p.2|) we introduced the convenient 
abbreviation (for Zk = Zk) 

H{u) = {u + Rk{u))-'^^Rk{u) 



(2 + |ln^.)^[exp(^j-lJ for (1.2) 



(3.3) 



{2 + fAnZk)Zke[u + Zkk 
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for (1.6) 



In the second, ghost-type, trace the function Hg{u) is defined similarly, but with a 
different factor Zk = 1. The last trace accounts for a possible difference between Zk 



and Zk in p.l|). Here H is obtained from H by replacing Zk by Zk and for r]_l we 



(2) 



may take equivalently Vx/ak or V^/ ak- 
in appendix A we display various trace identities for an evaluation of ( |3.2|) . 
They are exploited in appendix B to express the evolution equation for Wk as traces 
over simple operators. We note that the ultraviolet finiteness of the momentum 
integrals implied by the trace is guaranteed for the choice ( |1.2|) only if H{u) decays 
exponentially for large q^. This requires W'k > 0, Wk + B'^W'k > 0. Restricting the 
discussion to this case, one finds 
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f dxx2~^H{x) 
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with 



W+l 



■K^T 



(3.4) 



(3.5) 



The eigenvalues of the operator Vt are known explicitly . They are parametrized 
hy a {d — 2)-dimensional momentum g^ which "lives" in the space orthogonal to the 
1-2 plane, and a discrete quantum number n = 0, 1,2, ... which labels the Landau 
levels. The spectral sum for the function H{x) = H(WkX) reads 



Q-'Tt,,l[H{Vt)] 



d 



d-2. 



\d-2 



hi 27r to J (2^J 
■[{d-2)H{q^ + {2n + l)g\ui\B) 
+H{q^ + {2n + 'i)-g\ui\B) 
+H{q^ + i2n-l)g\ui\B)] 



(3.6) 



(The momentum integration is absent for d = 2.) Here 1^1,1 = 1, ..., N"^ — 1 are the 
eigenvalues of the matrix n^Tz in the adjoint representation. We note that for n = 



8 



and q^ sufficiently small the eigenvalue q^ — g\vi\B in the third term on the r.h.s. of 



[q) can become negative. This is the origin of the instability [14] of the Savvidy 



vacuum [|l| , which causes severe problems if one tries to compute the standard one- 



loop effective action in the background of a covariantly constant magnetic field |T5 
In our approach this problem is cured by the presence of an IR regulator. For d > 2 
eq. ( ^.61 ) can be rewritten as 



IT fri Jo 

■ Id f; H{x + {2n + l)g\ui\B) + H{x - g\ui\B) - H{x + gHB)] (3.7) 



n=0 



but it cannot be simplified any further in closed form. The other traces in ( p.4|) are 
given by 



iV2-l 



n-'TT^H{-D^)] = ^^ E 9Wi\B / dxx-^'^ ^ ^(x + {2n + l)gWi\B) (3.8) 

^ 1=1 -^0 n=0 

and similar for H and Hq- 

In this paper we use two different methods in order to (approximately) com- 
pute the spectral sums ( |3.7| ) and ( |3.8| ). In appendix E we shall represent them as 
Schwinger proper-time integrals [0. This method leads to compact integral rep- 
resentations which are valid for all values of B, but it has the disadvantage that it 
works only for a slightly simplified form of the cutoff function Rk{x) defined by (|1.6|) 
which leads, as we shall see, to ultraviolet problems. The second method consists of 
expanding the r.h.s. of (|3.7| ) in powers of B. It is applicable if gB <C fc^, and it works 
for any function Rk{x). The condition gB ^ k^ guarantees that we may express the 
sum over n by an Euler-McLaurin series, and that the ffist terms converge rapidly. 
In this manner (B.7|) turns into 



n-'TT,,,[H{W[{^-B')VT)] 



{N' - l)^mi^B')ri Hdx rdyxi-'H{x + y) (3.9) 

ZTT / JO JO 

Y.Ct[Y. ^r {aBf-mqB't-i j^ dx xi-^H^^-''\x) 



^ m=l \ 1=1 



with 



Ct = 7At(2''"-' - 1)^2™ - ,, ^ .,, (3.10) 

[2m)\ [2m — 1)! 

Here B2m are the Bernoulli numbers. In a second step one has to expand the B^- 
dependence of Wj.. Note that only even powers of B occur in this expansion. The 
group-theoretical factors S/Ii"^ ^f™ are discussed in appendix C. In particular for 
SU{2) these factors equal 2 for all values of m. 



If we use the Euler-McLaurin series (3.9) and a similar expansion for the trace 
J) in eq. (|3j), we find for SU{Ny. 

'"'-'-k'^ (iV^ - 1)4 - 2 5: r^E^4- l'^'^ I 



^ I - ;;^i \k' 



d 



with Trn defined in appendix C and the constants E^ given by 

i?™ = 7^(2'™-' - l)52rn. (3.12) 

The dimensionless integrals 

1 



'^ - ri Jo 



/ , \ 2m— 1 
d_2 a \ X 

dx x'2 -— 



\dx J e^ — 1 

I rl r'^^ rl TOO 'r2 

:z.fc2N-f / dxxi-^Hix) = dx 



2 — T] Jo Jo e"^ — 1 

r^^ = (ZkkY-' dx dyx-2-^H{x + y) 

2 — T] Jo Jo 

= fd, rd,^t!(£±j^ (3.13) 

JO Jo exp(x + y) — 1 

occur as a consequence of (|3.3| ) and the second equality uses (|1.2| ). For this choice 
we note for later use that in 4 dimensions 

ro''"^ = 52m-2, rt = 4 = 2a3) (3.14) 

where C denotes the Riemann zeta function. In contrast, the quantities rf and r| 
are not well defined for the choice ( |1.6| ). The ultraviolet divergence indicates an 
incomplete "thinning out" of the high momentum degrees of freedom for a simple 
mass like infrared cutoff. Even though eq. ( |3.11| ) was derived by choosing a specific 
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background field, this evolution equation does not depend on the background we 
used for the calculation. By employing derivative expansion techniques ||17|| , 



18 it 



should also be possible to derive (|3.11|) without ever specifying a background. In 
the case at hand the method presented here is by far simpler, however. 
In eq. ( p.ll ) we have introduced the anomalous dimensions 



V 



d 
"dt 



InZu 



V 



In Zu 

dt 



(3.15) 



A convenient choice for the wave function renormalization constants used in the 
infrared cutoff Rk is (fc„p denotes the transition to the nonperturbative regime) 



Vr^(O) for k>k 



np 



r\,n 



(0) for k<k 



np 



Zk — — 



1 
Oik 



(3.16) 



Here we account for the possibility that Vr^(O) may turn negative for k smaller 
than a "confinement scale" koo-, whereas Zk must always be strictly positive. More 
precisely, if |(ilnPF^(0)/dt| grows larger than one for small scales fc, we choose fc„p 
to be the scale where rj{knp) = 1- For k < knp all couplings run fast anyhow, and an 
improvement of the scaling properties of Rk by the introduction of a /c-dependent 
wave function renormalization seems not necessary. The choice Zk = a^^ guarantees 
that the infrared cutoff acts on the longitudinal modes in the same way as on the 
transversal modes. It implies that the flow equation for ^W^fc becomes independent 
of Ok in our truncation (except for an irrelevant constant)! In this paper we can 
neglect the running of a^. This can be infered from a first-order approximation to 
the solution of the general identities which govern the dependence of rA:[y4, A] on 
the background field ^[1]. On therefore has fj = 0. 

It is also convenient to express the flow equation in terms of renormalized di- 
mensionless quantities 



9^ 

Wk{^) 



d-Ary-l 



Zk'f 



g^k~'^Zk9 
g'k~^Wk{6) 



(3.17) 



Switching to a notation where dots denote derivatives with respect to "^ instead of 
6 and d/dt is now taken at fixed -^ one obtains 



d_ 

dt 



Wk{^) = -(4 - r])wk{^) + Mwki^) 
,((d-l)id-2) 



+ (2 - v)v,9'{wkm-^ ' T ' iN' - ly. 
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^'''^'''^^^ -{d-2)± rUCi - Ejrn^Mm' 



oo 

+2(rf - 2)vdg^ Y. r^E^4''^{2dr + const (3.18) 



m=l 



This is a nonlinear partial differential equation for the function Wk{'d) depending on 
k and d. One needs in addition the running of the renormalized gauge coupling g 
and r/, which are related by 



P^, = -i- = {d-A + r^)g' (3.19) 

For k > knp we have by definition Wk{0) = 1 and r] can be determined by 
^^,(0) = = r]-2id- 2)v,4'\,g' ((2 - r])Ct - (4 - r])E,) 
-(2 - ry) W (24 + ^(^ - 1)(^ - 2) (^2 _ ^yA ^^(0) (3_2o) 



With n = N,Cf = ^-2,Ei = ^ this yields 



jVd{d - 2) (26 - d)4'g^ - 2/1,(72^2 



"6 



Vdid - 2) (25 - d)4'9' + hag^w^ (3.21) 



with 



rf(rf-l)(rf-2),^. 



/i, = t;, I 2rf + ^ ^ '-{N' - l)4\ (3.22) 

W2 = ^fc(O) (3.23) 

For general d the constants Tq' , rf and rg depend on the precise choice of the infrared 
cutoff except for d = 4 where rg' = 1 is cutoff independent. The running of the 
renormalized gauge coupling g is now determined (cf. ( p.l9| )). It depends on the 
additional coupling W2 (|3.23|) which will be discussed in more detail in the next 



section]^. Specifying the initial value g'^{A) and the function w\{i}) at some high 
momentum scale A the form of Wk{^) and g'^{k) are completely determined by the 
flow equation ( |3.18D . Solving for A; ^ the function Wo{^) specifies the effective 



action in our truncation. If necessary, one has to replace for k < knp eq. ( p.21|) by 
1] = 0,g^{k < knp) = g^{knp). 

Before closing this section, we briefly comment on the range of convergence of 
the Euler-McLaurin series in our case. For large m one has 



lim Ci = d lim Em ~ vr-^™ (3.24) 



-2m 



4 



For W2 = we recover the result of ref. j^] except for the factor (25 — d) in the denominator 



in ( 3.21 ) which was (24 — d) previously. This difference is due to a slightly different choice of the 



Z-factors in the infrared cutoff. 
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For iV = 2 and rf = 4 we find (cf. i^T§) that the coefficients of (2^)'" in eq. (|XTB| ) 
diverge ~ '7r~^™i?2m-2 ~ 7r^^™2^^'™(2m — 4)!. For small nonvanishing -d the first 
terms of the series have a very rapid apparent convergence, but the series finally 
diverges due to the factorial growth ~ {2771 — 4)!. We can therefore safely use this 
series only for the derivatives w^'^\'d = 0) with finite n where convergence problems 
are absent since only a finite number of terms in the sum contributes. The situation 
is probably similar for N > 2 and/or d 7^ 4 as well as for many other choices of 
the infrared cutoff. In contrast, the original sums over n in eqs. ( |3.7| ), ( ^.8[ ) always 
converge since for B > the contributions from sufficiently high values of n are 
exponentially suppressed. An explicit evaluation of these sums is possible for the 
simpliffed masslike IR-cutoff ( |1.6| ). This is described in appendix E where we will 
also see the reason for the ultraviolet divergence of rf and r2 for this particular 
cutoff. 



4 The F^ approximation 

In order to get a first idea of the physical contents of the evolution equation ( |3.11|) 
we further simplify the problem in this section: We include in Wk{0) only terms 
which are at most quadratic in 6: 

Wkie) = Woik) + Wiik)0 + ^W2ik)e'', 

WAk) ^ l^-^J W,{0) (4.1) 

Thus the truncation for F^ is parametrized by three couplings: 

r,[A,A] = J d'^x {Wo{k) + ^W,{k)F^^Fr + ^W2{k) (F;,Fr)'} (4.2) 

Apart from the constant Wq, we keep here only the "marginal" coupling Wi and 
the lowest "irrelevant" coupling consistent with the truncation, W2- We concentrate 
on d = 4 and employ the infrared cutoff (1.2) with ( |3.13| ). The short distance or 
"classical" theory is specified for fc = A by {Z\ = 1) 

W^iA) = 1, Wo{A) = W2iA) = (4.3) 

In terms of the couplings Wq = w{{} = 0),g'^ and W2 ( |3.17| ), ( |3.23| ), one has for 

Woik) = k^w^{k)/g\k) 

W,{k) = f/g\k) 

W2{k) = tw2{k)/g\k)k' (4.4) 
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or 



g\h) = -g\ w7o(A) = w^2(A) = (4.5) 

with g'^{k) determined by eqs. ( p.l9| ), ( |3.21| ). For k < knp we use g^{k) = g'^{knp) 
and 



Wi{k) 
Wi{k) 



f 



g yf^np 



-Wi{k) 



(4.6) 



The flow equations for Wj{k) 



w 



ij) 



{d = 0) follow from appropriate partial differen- 



tiation of eq. (3.18) with respect to -(9, as, for example 

d 

-—w = riw + 4^m) 
ot 

{oo 
m=l 



d,m I 



•2\m-l 



+r 



w + 2'dw 



+-d{d-l){d-2){N^ 



{w + 2my 



d 



I / \2'Qw ) 
w{w + 2§w) 



{2mw'^ + (4m — d)'dww) 



w 



+A{d - 2)vag^ Yl rnTmE^rt''^{2d) 



\m—l 



(4.7) 



m=l 



We observe that wq does not appear on the r.h.s. of the evolution equations for 
Wjij > 1, and we omit in the following this irrelevant constant. For k > knp the 
evolution equation for the running gauge coupling reads {d = 4) 



dg^ 
dt 



247r2 



;ilA^ - 37^4^2) 



327r^ 



(7iV - 2H4W2) 



(4.8) 



whereas for k < knp one uses the flow equation for wi 



d 



g\k. 



npj 



<,TX^ 



. 3 'wf, 



(4.9) 

O). It 
until at 



where H^ = Wjc^h^ = rf + 3{N'^ - 1)4 = 2{3N^ - 2)C(3) for the choice 
is obvious that for W2 < the gauge coupling g'^{k) always increases ( |4.8| 
k = knp the anomalous dimension t] reaches one. If W2 remains negative for k < knp, 
the coupling wi decreases until it reaches zero at the conflnement scale fcoo > 0. (If 
we deflne g'^{k) = g^{knp)/wi, this coupling diverges at the conflnement scale.) The 
issue is different for W2 > 0: The coupling Wi would not reach zero since for small 
enough Wi the second term in eq. ( ^.9| ) would cancel the flrst term. One therefore 
needs an estimate of W2{k). 
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Let us first consider k > knp where the evolution equation for W2 reads (with 
p — _ 7 ^^4 _ _ 67 N 

-^2 - 720' ^2 — 180^ 

9 /. X 9^ ( 42/127 29 \ ,^ ^^, 

+ (2 - ,,) (5rJ + 5(iV2 - l)r5) u| - (2 - ,) (rj + ^-(N' - l)r«) t.3} 

We observe the appearance of the couphng W3. This is where the truncation of this 
section becomes relevant: We simply put w^ = here. We then observe that all 
remaining terms in the curly bracket in ( 4.10| ) are positive. Starting from W2(A) = 



we conclude that W2{k) is negative for all k between knp and A. In the perturbative 
region, where Ng'^/lGn'^ is small, it is easy to see that W2 is of the order g'^. In fact, 
we may neglect the terms ~ w| and ~ 77 in the curly bracket in ([4.10| ). One finds 
for the ratio ^ an infrared stable fixpoint 



9" 



dt\g^j~ g^ 3607r2"^' 



^hr =4-T + 7T77^^2ro' (4.11) 



-(*' - -li^-" v(*) - -i^ii^ ("2) 



which is approached very rapidly. (The last equation in ( [4.12| ) is for the choice (|1.2|) 
and N = 2.) It is interesting to insert the value (|4.12|) into the /3- function for g'^ 
( [4.8| ). Expanding in powers of g^ one has 

dg' 22N g' '^^iV^^^e ...... ^' 

-m = -^16^ - ^ (16^ + ('^ - 2)^(3)^^.2 

22 ^^iV /77 127 / 2 \\ ^^ ,^^„, 

= -yie;^ - ly + 1^^(')"^ I' - 3iv^j J (16^ ('-^'^ 

We note that without the term ~ W2 the coefficient ~ g^ exceeds the perturbative 
two- loop coefficient ~^7i^~T)7 by a little more than 10%. It is recomforting to 

find for the choice ([T^ ) [tq = |) the contribution from W2 in the same order of 
magnitude as this difference. This gives the hope that the dominant nonperturbative 
physics is already contained in the F^ approximation, whereas additional invariants 
(like F^ in this section) give only moderate modifications for k > knp. We emphasize 
that a full computation of f]g2 in order g^ should take additional invariants into 
account, as for example {FF^ or {D^F^'^Y. We also observe that for k = knp one 
has approximately Ng^/lGii^ = ^ such that the validity of perturbation theory 
extends roughly to all k > knp- 

The infrared fixpoint in W2/ g'^ implies that the coefficient W2 in ( |4.1| ) diverges 

W2r^-^ (4.14) 
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Because of the infrared divergence for k ^ 0, a result of this type could never have 
been found in standard perturbation theory. It could, however, be derived using 
the "-^ ~ ^"-approximation of the effective average action which we displayed in 
eq. ( |1.4| ). For this purpose one can neglect the Wj! terms on the r.h.s. of (|L 



and approximate Wl. = Wi. Then, with (2^) inserted into (|1.4|) , one obtains for the 
fc-dependent terms in F^ 



rk[A,A] = ^Tr,,iln 
-Tr^c In 



WiVt + {wi )D ^D + Rk{w{DT + (wi )D ® D) 

-D^ + Rki-D^)] (4.15) 



If one extracts the F^^^-term from these traces, one finds a (renormalized) coefficient 
which equals exactly 1^2, as extracted from the fixpoint (^4.12|) . Clearly ( |4.iq ) goes 
beyond the "-^ ~ ^ approximation. The terms proportional to r] and to iff could 
not have been obtained in this approximation. 

Let us finally consider k < knp. The evolution equation for W2 

d . , 9\knp) f 4,2/29 2 7 

:W2=4u;2 + ^-^ T-2ro' T^w^i 



dt 47r2 I ^ '^ V20 ' 180. 

+ (Srt + 1{N^ - l)A ^ - (rt + hN^ - l)r|) f^ \ (4.16) 



again contains only positive terms in the curly bracket if w^ is neglected and wi is 
larger than -i^. (For "U^i < 261 ^^^ positive term ~ w^/wi dominates largely.) The 
coupling W2-, which is negative at the scale knp-, will therefore remain negative for 

The investigation of this section strongly indicates a negative value of the cou- 
pling W2 for all scales. We conclude that the vanishing of Wi at a nonzero confine- 
ment scale fcoo is unavoidable within our truncation. This is a clear indication that 
the F"^ term in the effective average action changes its sign at the scale k^ol 

5 Gluon condensation 

If the coefficient Wi in front of the F^^F^'^ term changes sign for small enough /c, 
the solution A^ = Q cannot correspond anymore to the minimum of F^. Since F^ is 
bounded from below there must be a new solution with F^^F^^ > corresponding 
to the true minimum of F^ [Q. One therefore expects a phenomenon of "gluon 
condensation" , i.e. the minimum of Fq should occur for nonvanishing F^^F^^ . In our 
truncation this suggests that the minimum of Wk{9) should shift to a nonvanishing 
9q for small enough k. On the other hand, a negative value of W2 excludes a 
scenario where the minimum of Wk remains at = for all k > k^o, and then ^0 
continuously moves away from zero for k < k^o. A second local minimum of F^ 
must already appear for ki > k^ and become the absolute minimum for k < k^, 
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with koo < kc < k. At the scale fcoo the local ininiinuni at ^ = finally disappears. 
This scenario is closely analogous to a first order phase transition as a function of 
temperature. (In fact, the jump of the minimum of r^ as a function of k corresponds 
most probably to a jump of Fq as a function of temperature and therefore to a first 
order phase transition of the high temperature Yang Mills theory.) In order to 
describe such a "first order scenario" one has to keep at least the term ~ W^O^ 
in a polynomial expansion of Wk- We will show in this section that a polynomial 

approximation 

1 1 

Wk{^) = wii3 + -W2i3^ + -w^id^ (5.1) 

2 

is indeed sufficient for a qualitative description of gluon condensation. 

The first necessary ingredient is a positive coefficient w^lk). With (|1.2| ) we find 
in four dimensions the following /5-function for k > k^p (for details see appendix F) 

dt ^ " 167r2 1 30 V315 210 7 ^ ' 

+^(2 - r])r2W2 - 6C(3)(2 - vW2N^ + 23)wl - (m^ + 8)^2^3 + N^wA 

We neglect the term ~ w^ consistent with our approximation. For small values of 
g^ there is an infrared stable fixpoint of the ratio w^/g'^ corresponding to positive 

221r,g\k) 
"-^^) = 37800T6^ ^'-'^ 

Actually, all ratios Wn/g"^ reach perturbative fixpoints for n > 3 (cf. appendix 
F). This justifies the approximation i f).l\j at least for small enough g'^{k). (We 
observe Wn ~ g'^i'^~^)k~^^^~^\) A negative value of W3{k) would require that the 
positive term ~ W2 dominates the ws-independent part of (|5.2| ). This needs w^ > 
i2Af^+i9 60^3) ' or, for iV = 2, W2 ^ |. Such a large value of W2 does not occur for 
k > knp and we conclude that w^ is positive for this range of scales. 

For positive w^ the polynomial ( ^.1]) is bounded from below for ?9 > 0. For 
Wi > there is always a local minimum at -(9 = 0. Two additional extrema are 
present if 

wl > 2wiW3 (5.4) 

There is a minimum for positive -^o 



-W2 + Jw^ - 2WiWz 

Vo = (5.5) 

W3 



and, for u;i > 0, a maximum at 



~W2 - ^Wl - 2WiW3 

Vmax = (5.6) 

W3 
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For wi < the origin ^ = turns to a local maximum. The critical set of couplings 
where the two minima are of equal height {wk{^o) = 0) corresponds to 



3w^ 



wi = ^^ (5.7) 

For W3 smaller than the critical value, the absolute minimum occurs at '^q. This 
corresponds to the phenomenon of gluon condensation, with Euclidean action in the 
ground state given by Wk{^o). 

For k > knp where Wi = 1 a quick inspection of the flow equations (cf. ([4. 121), 



(|5.3|) ) shows that 2^3 remains larger than w|. There is therefore only one minimum 
at "i? = 0. For k < knp, however, wi decreases towards zero. In the region where w^ 
remains positive there is necessarily a scale ki where ( |5.4| ) is met and subsequently 
a scale kc where Wi(fcc) obeys ( ^.7|) . For k < kc gluon condensation sets in. Since 
kc > koo the establishment of the phenomenon of gluon condensation does not need 
the validity of the flow equations down to the scale fcoo- It is sufficient if the flow 
equations are a reasonable approximation for k > kc. 

What remains to be discussed in our truncation is the positivity of w^. The 
relevant flow equation for k < knp reads 

n /137 6 31 2\ 87 3 
30 IT05^^ + 315^0 + I^"^^^^^ 

^ + (12iV^ + 23)41| (5-8) 




w 



1 

Neglecting again W4 and for W2 < the only positive term in the curly bracket 
in (|5.8| ) is the one ~ w|/i(7p For k — *■ koo this term will flnally dominate the 
flow equation ( |5.8| ), drive w^ towards zero, and invalidate the truncation ( |5.1| ). We 
emphasize, however, that the critical condition ( |5.7| ) is always met at a scale kc 
where w^{k^ is positive. Also, from 

at IGTT^wf L V30 270 

,2 
2 



+ ?M[(9iv2 + Xy\ - {3N'' - 1)^] } (5.9) 




we learn that the term ~ W3, which was neglected in the previous section, is actually 
substantially smaller than the term ~ W2 for k between ki and kc- We therefore see 
a clear indication for gluon condensation within the F^ truncation employed in this 
section! 

The system of the three flow equations ( [4.9|) , ( ^.9|) , and ( ^.8|) (with W4 = 0) 



could be solved numerically with initial conditions given at the scale knp- It would 
be interesting to see if the expectation value of 9, i.e. 

Ooik) = Z^lk'g-'{knp)Mk) (5.10) 



has a tendency to stabilize at a fixed value before wi{k) or w^lk) reach zero. The 
corresponding value could then be taken as a rough estimate for the size of the gluon 
condensate. The breakdown of the approximations of the present paper for k — > k^o 
is, however, unavoidable for very general reasons that will be discussed in the final 
section. This excludes in the present approach a quantitative determination of the 
gluon condensate. 

6 Conclusions 

In this paper we have approximated the exact nonperturbative evolution equation 
for Yang-Mills theories [^ by a truncation where F^ is only given as a function 
Wk oi 6 = jF^j^yF^'^ . This ansatz is general enough to allow for a description of 
gluon condensation. A nonvanishing gluon condensate corresponds to an absolute 
minimum of Wk{0) for ^q > in the limit where the infrared cutoff k vanishes. 
We have not solved explicitly for A; ^ 0, but we have found strong indications for 
gluon condensation appearing already for nonvanishing k: First, the term linear in 
^ in a expansion of Wk{0) around ^ = turns negative aX k = k^o (sect. 4). This 
implies that for k < koo the "trivial" configuration F^j,^ = cannot correspond to 
the minimum of the Euchdean action F^. Second, we have found (sect. 5) that a 
new minimum of Wk{0) appears for k = ki, and that this minimum becomes the 
absolute minimum of Wk at the scale kc, with ki > kc > koo- The truncated set of 
evolution equations is qualitatively reasonable for all k down to kc and we interpret 
our results as a strong indication for gluon condensation. 

On the other hand, we note that for k in the vicinity of the confinement scale k^o 
our flow equations become very problematic. These problems are not originating 
in the general method but they occur rather as a particularity of our truncation. 
There are three different sorts of problems, and it is instructive to discuss them in 
detail: 

i) The flow equations are not valid for negative Wl.{9) = dWk/dO. For u = 
Wl.{9)x the momentum integrals implied by the traces of sect. 3, i.e. / dxxH(Wl^x), 
are ultraviolet finite only for Wl. > 0. In consequence, various terms in the flow 
equations ( [4.9| ), ( p.9| ), ( p.8| ) diverge for wi —* 0. The reason for this disease is the 
negative inverse propagator F^. ~ W^.q'^ for Wl. < 0. This implies in particular 
that the high momentum modes with g^ 3> k^ are all unstable. Even our optimized 
infrared cutoff ( |1.2| ) cannot remove this problem: The function H{u) defined in (O) 
grows ~ —u/k"^ for large {—u/k"^). 

The problem described here should, however, be viewed as a pure artefact of 
an insufficient truncation. The physics of the high momentum modes in four- 
dimensional Yang-Mills theory should be governed by asymptotic freedom, and F^ 
should be described by a term ~ ^Z^F^^F^'^ with positive Z^ for all k as far as 
the high momentum modes are concerned. In our truncation we have neglected the 
expected momentum dependence of Z^. A more reasonable approximation for the 
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term quadratic in the gauge field A^ would be of the sort 

\jd'xF^^Z,i~D'[A])Fr (6.1) 

with Zfc depending on the covariant Laplacian in the adjoint representation. For the 
functional form of Zklq"^) we expect for large q^ a /c-independent positive function 
Z{q^). (In lowest order the logarithmic dependence of Z on q^ should be deter- 
mined by the one-loop /9-function for the running gauge coupling. Indeed, external 
momenta of the gluons act as an independent infrared cutoff. For k"^ <^ q^ the 
running of Zk{q^) with k should essentially stop whereas it is given by the one- loop 
/3-function for k"^ ^q^-) In contrast, our truncation identifies Zk{q^) with a constant 
Zk{q^ = 0) = W^iO = 0). This leads to unphysical problems for negative Wj.. In 
reality, the quantity Zk{0) may well turn negative without affecting the high mo- 
mentum behaviour of Z}.{q^). We conclude that with a better truncation the high 
momentum modes are always stable if the momentum dependence of the propagators 
is properly taken into account. This implies positive u for x — *> oo and no ultraviolet 
problem can appear for the exponential cutoff ( |1.2| ). The present truncation can 
nevertheless be employed for all scales k where the inverse propagator for q^ th k^ 
can be approximated by Zk{0)q'^. One may hope that Zki^q"^ = k^) ^ Wl,{9 = 0) 
remains valid in the vicinity of kc- 

ii) For k < kc the absolute minimum of Wk occurs for ^o > 0. For this range of 
scales it would seem reasonable to expand VF^ around 6o rather than around 0. A 
polynomial expansion around the true minimum has turned out to be an important 
improvement for all models investigated so far with the help to the average action. 
The reason is that the relevant mass terms in the "ground state" appear then directly 
in the flow equations. Instead, an expansion around ^ = involves the "unphysical 
mass terms" corresponding to a metastable configuration. At the minimum of Wk 
one always has 

w;,{eo{k)) = (6.2) 

WHiOoik)) = A > (6.3) 

Differentiating (|6.2D with respect to k yields the running of ^0(^)5 i-e. 



This equation, together with flow equations for the couplings W^^-^Oq), could replace 
the flow equations of the last section for k < kc- 

Unfortunately, a quick inspection of eq. ( p.7|) reveals that the r.h.s. of such flow 
equations would diverge in the present truncation. It is instructive to investigate 
more in detail the contributions from the different modes to the evolution equation 
( |3.4|) : First we have modes with a positive "transverse kinetic term" (VF^(6'o) + 
B'^Wl!{9o))x = XB^x that we may call of type I. (Here x denotes the square of 
the transverse momentum.) The contributions from type I modes as well as from 
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ghosts are unprobleinatic for -t} = "^o- The type II modes have an inverse propagator 
Wl.{6)x + c with c > 0. Contributions from these modes suffer from the "high 
momentum disease" since WI.{9q) = 0. The corresponding ultraviolet divergences 
in the flow equations could be cured by employing Zk{q^) as discussed above. The 
infrared behaviour remains without problems. 

Finally, there could be modes of type III with inverse propagator VF^(6')x — c, c > 
0. This would lead to an infrared instability for the modes with x ^ due to a 
negative masse like term — c. Such an infrared problem cannot be removed by a 
modified truncation for the high momentum modes. It would be an indication that 
even at 9q the particular configuration with constant magnetic field introduced in 
sect. 2 does not correspond to the minimum of r^. This is obviously not possible 
within our truncation and we find correspondingly that no mode of type III exists 
in the spectrum for 6 = 6q. An expansion around "i^o becomes therefore possible for 
an enlarged truncation of the type (|6.1|) . 

We should mention that F^ is, of course, not only defined for the ground state 
configuration. Negative eigenvalues of F^ for a configuration like ( p.7|) are harmless 
as long as fc^ is much larger than the absolute value of a negative eigenvalue. Only 
once /c^ becomes of equal size strong renormalization effects set in. They lead to a 
"flattening" of F^. This is related to the general property that F^ becomes a convex 
quantity for k ^ 0. We are not dealing with this "approach to convexity" |T^ in 
the present investigation. 

iii) We observe that the spectrum of small fluctuations around the constant 
magnetic fleld conflguration ( p.7|) lacks Euclidean 5*0(4) rotation symmetry. This 
is not surprising since F^^ singles out two of the space directions. In addition, the 
spectrum is partially discrete - continuity exists only with respect to the transversal 
momentum. The lack of full rotation symmetry and the partial discreteness of the 
spectrum are actually related: For an 5*0(4) symmetric spectrum the continuity in 
two momentum directions must extend to all momentum directions. One expects 

(2) 

then a spectrum with a few separate particles, each of them having F^ depending 
continuously on a generalized 5*0(4) invariant of the type q^q^. The true ground 
state should be invariant under a generalized version of Poincare transformations 
[^] and therefore have an 50(4) symmetric spectrum. In ref. we have proposed 
a candidate ground state ior d = 4 and gauge group SU{N), N > 4, but no 
realistic candidate for the gauge group SU{3) has been found up to now. It would 
be very interesting to perform an analysis similar to the one presented here for this 
groundstate candidate instead of the conflguration (|2.71 ). 

In summary, our flrst attempt to investigate the phenomenon of gluon conden- 
sation with the help of nonperturbative flow equations is encouraging. Our simple 
conflguration (|2.7|) and the simple truncation give quahtatively the expected be- 
haviour of Ffc: The minimum of F^ does not occur for F^^, = for sufficiently small 
k, thus indicating the phenomenon of gluon condensation. The present truncation 
is, however, insufficient to describe the true ground state. In consequence, a detailed 
quantitative analysis of the size of the gluon condensate seems premature and should 
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wait for an investigation involving a configuration which is more similar to the true 
ground state. 

Appendix A 

In this appendix we collect a few trace identities which are needed in order to derive 
the evolution equation for Wk- We choose the covariantly constant background (f2.7D, 
( |2.8|) . Then eq. ( p.l3|) implies for any function / 



TT,,L[P\\f{D„P\\,P^)] = TT,L[f{d„l,0)] (A.l) 

because Trc[P||] = n^Uz = 1. Writing P± = 1 — Py and exploiting S oc P\\ it is also 
easy to see that 

Tt,,l[P^{D„S)] = TT,,L[f{D„0)] - TT,L[f{d„0)]. (A.2) 



>(0) 



Since F^^, is a constant matrix, the operator s^j,^ of ( 2.14 ) commutes with P^ and 
Pj. and satisfies d^s^i, = 0. This fact can be used to show that 

Tr.z.[/(Pf,Pf;s)] 
= Tr,i[/(0,l;s)] + Tr,[/(l,0;0)]-Tr,[/(0,l;0)] (A.3) 

If one subtracts the same expression with s = one obtains 

Tr.,[/(Pi°\p(°);.)-/(Pr,4°);0)] 
= Tr,i[/(0,l;s)]-rfTr,[/(0,l;0)] (A.4) 

In the last step we used that Tr^j^ = dTix for an operator ~ 6fj_u- In the above 
identities the function / may also depend on further operators provided they com- 
mute with those displayed explicitly and do not introduce any additional colour or 
Lorentz index structures. 

For the evaluation of Ui in eq. ( |B.5D we need another important relation: 

Tr^cLiPLfiVT)] = Tr,,,[/(-D2)] (A.5) 

It follows from the fact that the operator {VT)f_iu = —D'^^^iv + '^igF^u, when restricted 
to the space of longitudinal modes (a^ = {PL)'^a^), has the same spectrum as —D^ 
acting on Lorentz scalars. The proof makes essential use of the identity 

Vt{D ®D) = {D® D)Vt = -{D ® D){D ® D) (A.6) 

which holds true whenever the gauge field contained in the covariant derivatives 
obeys D^F^^ = 0. 
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Appendix B 

In this appendix we use the trace identities derived in appendix A in order to simphfy 
the flow equation ( p. 21) . The first trace on the r.h.s. of eq. ( p.l|) can be simphfied 
by using the identities of appendix A. Inserting a factor of 1 = P\\ + P± leads to the 
decomposition 

Tr.,Ji7(rf[A^])l=rJi+Ti^ 



with 

Ti = Tt^cl 
= Tr^L 
where ( |A.1|) was used, and 

= Tra;cL 



P\\H{WIVt + [Wl - —]D ®D + W'^S) 
H{-Wld^ + [Wl - —]d ®d + ly^'s 

Ok 



P^H{WIVt + [Wl - —]D ®D + W'^S) 

ak 

H{WIVt + [Wl - —]D ® D) 
ak 

H{-wid^ + [wl - —]d ® d) 

ak 



-Tra;L 

where ( |A.2|) was exploited. Let us write 

T:i^,L[H{Tf[A,A])\ =Ui + U2 
with Ui the "nonabelian" trace 



1 



H{WIPt + [W'k ]D ® D) 

ak 



and U2 the sum of TJ and the second term of ( p.3|) : 



U. 



Tr^L 
-Tr^;!, 



Tl p(0) 



>(0)l 



H[-d'[WlP'^^ + -Pr]+W';s 
ak 



V ak J 



:b.i) 



(B.2) 



(B.3) 



(B.4) 



(B.5) 



(B.6) 



It is quite remarkable that if we now apply the identity (|A.4|) to f/2, the longitudinal 
contribution drops out completely and the result becomes independent of the gauge 
fixing parameter a^: 



U2 = Ti,L[H{-d^Wl + W^s)] - dTi.,[H{-d^Wl)] 



(B.7) 
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The operators entering ( [B.7| ) are easily diagonalized in a plane-wave basis. A stan- 
dard calculation yields, for Wl, >0,W!, + B^Wj! > 0, 



n-^Uo = 2v, 



WL + B^W" Wi \Wl 



Jtt I I Jtt 1 ' / dx x^~'H(x) 



(B.8) 



with Vd = [2'^+^7r'^/^r((i/2)] ^. (As always, the argument of Wk and its derivatives is 



understood to be ^B^. 



Next let us simplify the trace f/i by inserting a pair of projectors: 



Ui 



Tr 



xcL 



PtH{Vt[WIPt + —Pl]) 



+ Tr 



xcL 



1 



PlH{Vt[WIPt + —Pl]) 

ak 



1:i^,l[PtH{WIVt)] + Ti.,l[PlH{—Vt) 

ttfc 



The a-dependence of Ui is contained in 
AUi{ak) = TtxcL 
= Tr 



(B.9) 



1 



Pl{h{—Vt) - h{w;,Vt)} 

ttfc 

H{-—D^)-H{-WID^ 

Oik 



(B.IO) 



In the last line of (|B.10|) we made use of the identity ([A.5|) . 

By a similar combination of the trace identities we can also evaluate the last 
term on the r.h.s. of ( |3.2| ) 

""M \ ak V "fc / J 









:b.ii^ 



Evaluating the second term in a plane wave basis yields 






-f"^^*"!^©-^!!;)) (^■") 



At this point we have exploited the various trace identities as much as possible. 
Combining these results yields the flow equation (|3.4|). 



Appendix C 

In this appendix we discuss the group theoretical factors J2e^e"^- The LHS of the 
evolution equation is dtWk- The argument of Wk is \F^^F^^ = ^B^ which is mani- 
festly independent of the unit vector n^ which specifies the direction of the field in 
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"color space" . The r.h.s. of the evolution equation consists of expansions such as 
(|3.9|) which involve the factors J2e^e"^- As {z/^} are the eigenvalues of n^T^, we can 



rewrite them as 



V- ^2m ^ ^.,^.. . . . ^^.^TrjT.^T,, ■ ■ ■ T,,J (C.l) 



where the trace is in the adjoint representation. The question is whether the in- 
variants (|C.1|) are all independent of the direction of n^. In appendix D we explain 
in detail that generically this is not the case. If the orbit space of the gauge group 
in the adjoint representation is nontrivial, different n's can lead to different sums 
J2e^£"^- The resolution to this puzzle is as follows. For m = 1 we can use the 
standard orthogonality relation 

Tre[T,T,] = NSy, (C.2) 

to prove that J2e ^e = Nn^Uz = N is independent of the direction of n. Likewise, if 
the symmetric invariant tensor Trc T(^zi ■ ■ ■ ^Z2m) is proportional to the trivial one, 
'^(2122^2324 ■ ■ ■'^22m-i22m)) ^^ cau agalu use the normalization condition n^riz = 1 to 
show that the r.h.s. of (|C.1| ) is independent oin. The situation changes if there exists 



a totally symmetric invariant tensor ^i22--22m which is different from the trivial one. 
Then we might have 



Tr, 



-'(21 ■ ■ ■ -'22m) 



TmO(2i22 ■ ■ ■ ^22m-l22m) "•" -'2i22---22m I'^^-'^j 



with some coefficient r^. (If there exists more than one T an appropriate sum is 
implied.) Ingeneraln^^n^^ 7^^^^ -y^^iUl^gfjij-gc^ion dependent [0. If some invariant 
tensor T exists, the correct way of deriving the evolution equation for Wk is to 
compare coefficients of a fixed tensor structure on both sides of the equation. Clearly 
the l.h.s., dfWkijF^^) , gives rise to the trivial tensor structure only. Therefore only 
the Tm-piece of ( |C.3P should be kept in (|C.1|) and the (n-dependent) part coming 



from T has to be discarded. Thus ( |3.9|) may be used on the r.h.s. of the equation for 
Wk provided we interpret J2 ^i"^ as the coefficient r^. On the other hand, nontrivial 
T permits to construct additional invariants from an even number of covariantly 
conserved Ff^^^. The truncation W{6) is not sufficient anymore to parametrize the 
most general effective action for constant magnetic fields of the type introduced 
in sect. 2 (with covariantly constant -F^ijy). The evolution equation for the new 
invariants can now be extracted by projecting the r.h.s. on the appropriate tensor 
structure. We will not pursue this generalization in the present paper. 

In appendix D we show that for SU{2) this complication is absent. There exists 
no additional invariant tensor T, and one finds the n-independent result 

E^r = 2 (C.4) 

1=1 

for all 171 = 1, 2.... 
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Appendix D 



2m 



In this appendix we investigate in more detail the group-theoretical quantities J2i ^i 
which occur in many calculations involving covariantly constant backgrounds of the 
type Afi^ = n^A^. Here we consider an arbitrary (semi-simple, compact) gauge group 
G with structure constants f^y^. For a fixed unit vector n^ we consider the matrix 

f^yz ^ n'"{T'^y' = iP'^'n'" (D.l) 

The numbers ui,! = l,...,dim G are the eigenvalues of n: h^^tpf = uiipf . This 
equation can be rewritten in a more suggestive form. Let t^ denote the generators 
of G in an arbitrary representation: [t^,ty] = if^y^t^. If we define 

n = nH\ tPi = ijft' (D.2) 

the eigenvalue equation becomes 

[n,'ipi\ = vi-ipi (D.3) 

Clearly the z/^'s do not depend on the representation chosen. We would like to know 
how the spectrum {z/^} depends on the vector n. First of all, it is clear that if V 
is any group element in the t-representation, the matrices ft and n' = VnV~^ have 
the same spectrum, i.e. the spectrum is constant along the orbit of G in the adjoint 
representation. If two directions n and n' are not related by a group transformation, 
then the spectra {//^(n)} and {//^(n')} can be different. Typically, for G large enough 
[ Pm , the orbit space is indeed nontrivial, and the spectrum "feels" the direction of 

Let us go over from the basis {T^} to the Cartan-Weyl basis {Hi, E^} of the 
abstract Lie algebra. Here a G R'^ are the root vectors and i = 1, ...,r = rank G. 
For definiteness we assume that the t^'s are in the fundamental representation where 
we write {hi, e^} for the Cartan-Weyl basis. Thus 

[Hi,E^] = aiE^ and [/ij, e^j] = ajC^ (D.4) 

We assume that the generators hi of the Cartan subalgebra are given by diagonal 
matrices. By an appropriate transformation n -^ VnV~^ any n can be brought to 
diagonal form. Therefore, in order to investigate the ra- dependence of {z//(n)}, it is 
sufficient to consider fi's which are in the Cartan subalgebra: ft = Yh=i ^i^i- For 
this choice 




[n, h,] = (D.5) 

and the nonvanishing eigenvalues ui = u^ are given by 

r 

i^d = ^niai (D.6) 



i=l 
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Therefore the quantities Y,i ^i can be computed exphcitly from the root system: 

2m 

I (D-7) 



iimG 




/rank G 


>>- = 


- >; 


2^ UiUi 


1=1 


Toots{a} 


\ i=i 



Let us consider a few simple examples. For G = SU{2) we have a = ±1. Thus 

^z.f- = 2, m= 1,2,3,... (D.8) 

1=1 

depends neither on the direction n^ nor on the power m. This degeneracy can be 
understood by noting that for SU{2) the square of the matrix ( p.l|) is the projector 
P±: hh = P±. This means that E; ^'f™ = Tr(Pf) = Tr(P^) = 2, as it should be. 
Contracting eq. ( |C.3| ) with n^^...n^^"' and comparing the result to ( |D.^ ) we see that 
there exists no nontrivial invariant tensor T and that r^ = 2 for all m. 

For G = SU{3) we have r = 2, and a 2-component unit vector (ni,n2) specifies 
the direction of the field in the Cartan subalgebra. Using the explicit form of the 
roots it is straightforward to derive that 

^ ^2„^ ^ 2^-2™ J^^^ ^ 73^2)'™ + (ni - VSns)'"^ + (2ni)2-] (D.9) 

For m = 1 and m = 2 it turns out that this expression depends on rii and n2 only 
via nl + n^ = I, and one obtains the direction-independent results 

E^f = 3, E^f = ^ (D.IO) 

1=1 1=1 ^ 

Starting from m = 3, the invariants are explicitly n-dependent. Writing rii = 
cos 9, n2 = sin 6 we find for m = 3 

* 3 r 1 

Y^yf = — ^11 cos^ ^ + 15 cos^ e sin^ ^ + 45 cos^ 9 sin^ ^ + 9 sin^ ^J (D.ll) 

As discussed in section 3, the n-dependence is related to the existence of a nontrivial 
invariant tensor '/^^..^g. However, we are not going to calculate the corresponding 
coefficient t^ here. 

Appendix E 

In the regime where gB /k^ ^ 1 the use of the Euler-McLaurin expansion ( p.9|) be- 
comes questionable and we should look for an alternative representation of the 
spectral sums ( |3.7|) and (|3.8| ). In this section we use the Schwinger proper-time 
representation |]16[. It can be easily applied only for the simplified cutoff function 



Rk{x) = Zkk\ (E.l) 
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^W-^^^f^'^Vo""--'""^-" (E-^) 



In this case one may write {x = T>t) 

dtRkix) d 

Inserting this representation into ( |3.7| ) we may employ 

OD -1 

Y: exp[-sW;,g\i^i\Bi2n + 1)] = _ . ,.,.,,, , „. 
^0 2smh[sWf,g\ui\B] 

and 



Jo 
as long as Wj^ > 0. One finds 



vj\2'-'n'-Hsw;,y- 



rr\l- 



C , „ , o^ -.-^ , , „ / as (4_d)/2 -sZfefc2 



.^(Z.F) g SI.IB^ ^. 



exp(-sVr^^|z/i|5) + exp{+sW;,g\ui\B) 



2smh{sWl^g\ui\B 



(E.3) 



where the last exponential is due to the unstable mode. For Z^^k"^ < Wlg\vi\B 
it makes the s-integration divergent at the upper (i.e. IR) limit. In conventional 
calculations of the one- loop effective action this creates a problem from the outset, 
because one attempts to work at fc^ — * there. In the present formulation everything 
is well defined for k"^ sufficiently large, and one interesting question is how the 
renormalization group fiow behaves as one approaches k"^ k. 9 from above. 

In the UV limit s — >^ the terms inside the square bracket in (|E.3| ) behave as 

' ' "^ - + Ois) (E.4) 



' ' 2g\vi\BWl{B'^/2) s 

While the 0(s)-terms do not lead to UV-divergences for d < Q, the term ~ 1/s leads 
to a divergent contribution to the proper-time integral. Though the factor g\i'i\B 
cancels against a similar one coming from the density of states, this divergent piece 
is still field-dependent because of the -B-dependence of W^. This UV divergence 
shows a failure of the truncation for the mass- type cutoff function Rk = Z^k'^. The 
latter may be used only together with the approximation W^l = on the r.h.s. of 
the flow equation. The divergent piece in the proper-time integral is an irrelevant 
constant then. 

Using ([E.3|) and a similar formula for the scalar traces in eq. ( p.4|) , we obtain 
(up to an irrelevant constant and for Z^ = l/a^jfj = 0) 

^W,(-B') = (47r)-^/2 V Q\u,\Bk^ r~ds s^-"'^ ■ 



dt '2 



(47r)-rf/2 y^ g\ui\Bk^ ds 
1=1 -^0 
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exp{—sk^) 
sinh{sg\iyi\B) 



LY- 



d-i 



2smh{sWlg\pi\B) 



+ 2s\Yi\i{sWlg\pi\B) 



-sZuk'^ 



(E.5) 



This evolution equation is the analogue of ( |3.11 ) with the additional assumption 



W' 



0. Contrary to the Euler-McLaurin series it is valid even for strong fields 



gB ^ P. 

For gB -C k"^ the r.h.s. of ([E.5|) can be expanded in powers of B. Apart from the 
different form of Rk, this reproduces the Euler-McLaurin expansion. Expanding up 
to order B"^ we find, for instance, 



4feU4^ + 



dt\g' 



9^ 



127 4,2 



(E.6) 



with Vq = 2. This result is the counterpart of eq. ( 4.11 ) which had been obtained 



with the exponential cutoff for which rg' = 1/6. In accordance with the ( p.l3|) we 



find an additional factor of 12 in the second term on the r.h.s. of (|E.6| ). Hence also 
the value of the fixpoint W2*{k) is 12 times larger than the result ( [4.12| ). In view of 
the discussion following ([4.13|) this means that for the mass-type cutoff Rk = Zkk"^ 
this higher-order correction is much larger than for the exponentially decreasing R^ 
of ( p..2|) . This is probably closely related to the ultraviolet problems and indicates 
that, though computationally more difficult to handle, the cutoff ( |1.2D should be 
used for reliable estimates. 



Appendix F 

In this appendix we derive the flow equation in the F^ truncation. We start from 
the evolution equation for w{'d) which follows from differentiating ( |4.7| ) with respect 
to?9 
d 



dt 



w = (4 + iq)w + Adw^''^^ 



-(2 - r^)v,g^w-^ (^ - 2) E ^m{Ct - E^)4 



d,m 



m=l 



(f 



{2^w^)"'~' ((8 - d)m - d)ww + ((4 - 2d)m -d+ —)^w^ + (4m - d)^ww 

+2(m - l){2^w^)"^-^w^{w + 2'dw){2mw + (4m - d)-dw) 

4^i;(3) + 2^«;(^) 4(w + ^'«;(3))(3m) + 2'dw^^^) + 2'dw{bw^^'^ + 2'dw^^^^] 



,(3) 



+r 



w + 2'dw 



[w + 2dwy 



{w + 2§wf 
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2v? + Mibw^^^ 3 , 
-a — -^- h -d- 



•dib 



1 1 t"J 



w 



{w^ldib) 2 u;2(^w + 2t9tl;) 



+ 2(i 



t9w2(3w + 2^ti;(3)) 
w(zi; + 2'dwy 



+-d{d-l){d~2){N'-l)ri 



w 



(3) rf + 2 



m)^ 



W7 



2 w;2 



+8(d-2)WE 



mlm 



l)r^E^ro''™(2^) 



rf,m/r) Q\m— 2 



(F.i; 



m=l 



Taking one further -(^-derivative at "i? = we obtain for rf = 4 and k > knp the flow 



equation for w^ (with C2 — E^ 



29 Ai4 771 



137 
10080 



and Er, 



31 
30240 - 



d 



9 



+ r])w3 
2 V442 137 



137 \ 43 87 42 

— r^j Tgro' + —{2 - r^JTsro' W2 



167r2 IV315 210 

(2 - ri)r^{3w4 - 5IW2W3 + 105wl) 

■3(2 - ri){N^ - l)4iwi - 9w2W3 + I2wl] 



{F.2] 



For the infrared cutoff (|1.2| ) one has Vq = |,ro' = 3^ 



^ Tq'" = — 3^ and we observe that the 
perturbatively leading term ~ g^ is negative. This implies a positive perturbative 
fixpoint value for the ratio w^/g^. 

This discussion can easily be generalized for arbitrary w„. In lowest order in g"^ 
the flow equation ([F.l| ) simplifies considerably {k > k„ 



^npj 



d_ 
dt 



w 



4w + Mw^^'^ 



-8{d - 2)v,g^ E m{m - 1)t^4'"' (ct - 2E^) (2^) 

m=2 

This implies for n>2 the flow equations 

^Wn = e-^Wn -{d- 2)t;rf(722"+in!r„ro"'" (ct - 2E„) 
and the infrared fixpoint values 



m-2 



(F.3) 



\id-2)v,g^J^ 
For d = 4 and with ( p.l4| ) this yields 



23-"n!r„ro''" {C^ - 2E, 



[FA] 



(F.5) 



Wr, 



rf. 



9^ 



dr, 



n67r2 

24-n^l 



(2n-l)! 



TnB2n-2 



22n-i _ ;^ 

2n 



5 



2n 



(F.6) 
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where for SU{2) 

2d2 ^ _127 2d, ^ TZ]_ 

r2 270' ra 37800 ^ ' ' 

The series of dn is alternating as long as the bracket is dominated by -1. 
Finally, for k < knp the flow equation for w, follows from ( H*'.!! ) as 

(9 „ g^ ( 4 3/137 6 31 2\ 174 42 , 

-6rt (^ - 17^ + 3541 

\Wi w( WiJ 



4^ 

Wi Wi wf 



-6(iV=- 1)4 1^-9^ + 12^U (F-8) 
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